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Abstract. - The generalization of 4D confining string theory to the S'l7(3)-inspired case is 
derived. It describes string representation of the Wilson loop in the SU (3)-analogue of compact 
QED extended by the #-term. It is shown that although the obtained theory of confining strings 
differs from that of compact QED, their low-energy limits have the same functional form. This 
■ fact leads to the appearance of the string #-term in the low-energy limit of the SU (3)-inspired 

confining string theory. In particular, it is shown that in the extreme strong coupling regime, 
the crumpling of string world sheets could disappear owing to the string S-term at 9 = tt/12. 
' Finally, some characteristic features of the SU(N)-case are pointed out 

One of the ways to construct string representation of QCD is to use the method of Abelian 
projections |jj (see ref. Q] for recent reviews). The assumption of monopole condensation 
employed within this method has recently been proved with a high accuracy by the lattice 
calculations in ref. ||. To model the condensation of Abelian-projected monopoles one can 
treat them as a grand canonical ensemble with the Coulomb interaction. In the 3D SU(3)- 
inspired case Q , this led to a certain generalization of the theory of confining strings j^j , which 
was originally invented for the description of confinement of electric charges in 3D compact 
QED ||. The aim of the present letter is to extend the results of ref. || to the realistic 
4D case with the inclusion of 0-term. The model in which confining strings will be studied 
is thus nothing else, but the 5J7(3)-version of 4D compact QED with 6>-term. In this way, 
we shall see that similarly to the case of compact QED 0, in the S'[/(3)-inspired model the 
field-theoretical 6*-term generates eventually string 0-term. At a certain critical value of 9, the 
latter one might be important for getting rid of crumpling of the string world sheet. 

Consider first the case without the 9-term. Within the so-called Abelian dominance hy- 
pothesis H , the S'?7(3)-inspired theory under study describes free diagonal gluons and Abelian- 
projected monopoles (c/. the 3D-case Q,|j|). The action of this theory can be written as 
follows 



( x ) Throughout the present letter, all the investigations will be performed in the Euclidean space-time. 
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\ J d\¥l v + \J d 4 xj d*y}™»(x)D (x - y)j™»(y). (1) 

Here, F M „ = <9 M A„ — d„A^ is the field strength tensor of diagonal gluons = (A^, A^) and 
Dq(x) — 1/(Att 2 x 2 ) is the 4D massless propagator. Besides that, in eq. (Q) the iV-monopole 
current j™ on (x) is defined for N > as follows: 

N 

j™ n (x)=g m J2<lc, a fdzZ(T)5(x-x a (T)), (2) 

a=l J 

whereas for N — 0, j™ on (a:) = 0. Here, the magnetic coupling constant g m is related to the 

QCD coupling constant g by the equation gg m — Air. Next, in eq. (||), we have parametrized 

l 

the trajectory of the a-th monopole by the vector iJJ(t) = j/° + z^(t), where y" = J ' cZrx™(r) 

o 

is the position of the trajectory, whereas the vector zJJ(r) describes its shape, both of which 
should be averaged over. Finally, the magnetic charges of monopoles are described by six root 
vectors of the group 5(7(3): qi = (1/2, V3/2), q 2 = (-1, 0), q 3 = (1/2, -y/3/2), q_ Q - q«. 

The 0-term by which we are going to extend the theory (|l|) originally having the form 
"iff 7 / d^xG^vG^v can be rewritten modulo full derivatives as J d*a;A / J? oa . Here, 
Ofiv — ^£fj.vXpOx p , G MI/ = F M „ + FJf° n , and the monopole field strength tensor, F™° n , is 
related to the current j™ on according to the modified Bianchi identities: c^F" 1 ™ = j™ on . 

Let us now concentrate ourselves on the derivation of an effective field theory describing the 
monopole ensemble. The desired monopole part of the partition function, Z nlon [A M ], (which 
should then be averaged w.r.t. the action | J d^xF 2 ^) is nothing else, but the statistical 
weight of the TV-monopole configuration, 



2[j7 n ,Aj =exp -\ J d 4 x J d*yj™»(x)D (x-y)3™»( y )- 1 -^ J d 4 xA,j 
summed up over the grand canonical ensemble of monopoles: 



mon 



(3) 



oo tN 



mon 

iV=l 



Here, ( a e Smon is the fugacity (Boltzmann factor) of a single monopole loop, which has the 
dimension (mass) 4 . As far as the action of a single a-th loop is concerned, it obeys the estimate 



•Smon c< J dry , where the integral expressing the length of the loop is supposed to be 

of the same order of magnitude for all the loops. Next, the average over monopole loops in 
eq. (0) has the following form: 



A' 



(° [jri > = n / d v / vz a n \z a \ Yl ° K 1011 ] (5) 

a=]V J a a =±l,±2,±3 

In eq. (^) , /i [z a ] stands for a certain rotation- and translation invariant measure of integration 
over the shapes of monopole loops. Similarly to the case of compact QED flC] , the concrete 
form of this measure is irrelevant to the final result for the effective action. It is only im- 
portant that this measure is normalized in the standard way, (1) = 1, and from now on this 
normalization condition will be implied. 
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Equation 



can further be rewritten as follows: 



^ = E / v '^ z [jv. a mJ <<s Op - j" lon )) = 

N=l ' J 



^ A/! 

JV=1 



exp 



mon 



(6) 



The sum over AT here is equal to 



exp 



a=l 



exp 



£>z^[z] cos (g m q a <j) dz AI (r)l AI (x(r)) 



2c / AE cos 

0=1 



A 



(7) 



The derivation of this formula was based on a certain evaluation of the integral over shape s of 
monopole loops in the strong coupling limit j> 1 analogous to that of compact QED [|lc' , 1 1 



Besides that, it has been assumed that characteristic distances between loops y (y°) as well 
as characteristic loop sizes \J [zfy 2 are of the same order of magnitude for all a's. The 
respective UV momentum cutoff A (which is much larger than 1 / y (z®) 2 ) obeys the estimate 
A 



gy yyfy / (z°) . Finally, in cq. (Q), |Lj| denotes the modulus of the vector only w.r.t. the 

Lorentz indices, but not w.r.t. the Cartan ones, i.e. |1 M | = ^ 1 1 M | 1 , |1 M | 2 ^ = (yjtfjjii sftfjty ■ 

We should now carry out the integration over the Lagrange multiplier 1 M stemming from 
eqs. @ and (Q). Since the field 1 M has no kinetic term, this can be done in the saddle- 
point approximation. For every Cartan index a = 1,2, the respective saddle-point equation 



iA_ -a 
2C J V 



It can be solved w.r.t. q Q |1 M | by adapting for Z£ 



reads JL £ q a a sm(^M) 

' M ' Q = l V ' 

3 

the following Ansatz Z" = |Lj|° j"/ and representing |j M | a as £ 9aJ^' where j^ 1 ^ = 

(|j M r/V3+ Ij^l 2 ) / A J (2) = -2|j M | 1 A J (3) = (|j/A/3- UpI 2 ) /V3. The so-obtained 
solution reads q a |1 M | = —iA arcsinh (Aj( Q )/(2£)) ■ By making use of it and performing the 
rescaling of fields <?A M = A" ow , gj^ = j" ow we eventually get from eqs. (|TJ), (||), and (^) the 
following expression for the full partition function: Z = J I?A AI I?j AI e _ ' s , where 

s = 4? / dixF ^ + ^I dix I di y^w D °( x - yy^y) + ^ / ^-m* + ^ [W ■ ( 8 ) 



Here, 



V) arcsinh (^«) _ W 1+ 



(9) 



( 2 )Note that similarly to the case of 4D compact QED fij|ll|, cq. (fij) can also be derived by using the lattice 
regularization. 
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is the multivalued potential of monopole currents. 

As we shall see below, it is convenient to unify the kinetic term of the field and the 
interaction of monopole currents by introducing the following antisymmetric tensor field (else 

^ + h M „. Here, the field h M1/ is unambiguously 



called Kalb-Ramond field |12|) P) B M „ 



related to the current as h 



\ p df J d 4 yDo(x — y)j P {y) and thus obeys the modified 



Bianchi identities d^h^ = j„. ( ) In terms of the field B M „, the action (g) takes the following 
more compact form: 



— / d 4 xB 2 



32tt 2 



(1 ^B^jyB^jy 



V 



(10) 



In what follows, we shall be interested in the string representation of the Wilson loop 
defined as (W(C)) = 3 (trP exp ( i § dx^A^T ) V Here, A^ ot is the total vector potential 



which includes also the monopole contributions, and T = (2 '2) w ^h denoting the 
Gell-Mann matrices. Analogously to the case when monopoles are absent, Stokes theorem 
yields for the Wilson loop the following expression: 



(W(C)) = i(n-«-xp 



(11) 



Here, £^[2;, S] = J da^(x(£,))S(x — x(^)) is the vorticity tensor current defined at a certain 

S(C) 

surface S(C) bounded by the contour C and parametrized by the vector £p(£) with £ = (£*, £ 2 ) 
standing for the 2D coordinate. In a derivation of eq. ([ll]), we have omitted the path ordering, 
which is possible due to the fact that both A3 and As are diagonal. Note that by virtue of 
the quantization condition gg m — 4tt one can prove that eq. ([ll]) is indeed independent of the 
form of the surface S. 

The average in eq. (|ll]) is first performed over the free part of the A^-action, J d 4 xF^ ll/ , 
after which the result should be weighted with Z ^j™ on , an d summed up over the grand 

canonical ensemble of monopoles in the sense of eqs. (§)-(^)- This procedure can be simplified 
by rewriting the Wilson loop in terms of the above-introduced field B Mt/ as follows: 



(W(C)} 



exp 



d a;B„„Q a S] 



(12) 



Here, (. . .) stands for the average w.r.t. the action (|l0|), and the vectors Q Q 's, which denote 
the charges of a quark of the a's colour w.r.t. the diagonal gluons, have the following form: 

Ql = (~2> iTs) ' Q2 = (s' iTs) ' Qa = (°' ~7s) ■ Besides that > in the derivation of eq. © 

3 

we have used the identity tr exp (iaT) = ^2 exp (iaQ Q ) valid for an arbitrary vector a. 

In the low-energy limit, Ij^l <C sC/A, the monopole potential (B) becomes a quadratic 



functional. This yields the kinetic term of the field B M „, V 



Y2iy 



J^HLa: Where 



V = and 
limit takes the form 



d f _ l Ti 1 j\+d\Ti f _ l v + d u 'B\ f _ l . The Wilson loop ( |12| ) in such a low-energy 



( 3 ) In the formal language, this equation represents the Hodge decomposition theorem. 

( 4 ) Clearly, the same Legendre transformation, which made from the current jJf on the dynamical field j M , 
makes from F™° n the dynamical field h^. 
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(W(C))\ 



low cn. 



-'low cn. 



a=l 



^2 / eX P 



d 4 X 



it 2 



B„,,B,, 



4<? 2 32tt 2 ^ ^ 2 



(13) 

where i?i owcn . is given by the second line of eq. (|l^) with £ M „ set to zero. Carrying out the 
B^-integration, whose details are outlined in the Appendix A, we get the following result: 



Here, 



(W{C))\ 



low cn. 



exp 



12 



2g 2 <j> dx^ <j> dy^Dn^x - y)- 



c c 



d 4 yD m (x - y) ( E J[i „(a;)S /JI/ (y) + -^■'E IXI/ (x)E flv (y) 



(14) 




(15) 



is the mass of the field B M „, and D 



(x) — mKi{m\x\) / (47r 2 |x|) is the massive boson propa- 
gator with K\ standing for the modified Bcssel function. 

Upon the derivative expansion of eq. ( |l4| ) (analogous to the expansion performed in ref. Jl3| 
within the stochastic vacuum model of QCD), one gets as a few first local string terms the 
usual Nambu-Goto and rigidity terms, responsible for confinement of electric charges and 
stability of strings. Besides that, this expansion yields string 0-term [|l4| [^)] equal to icv, 
where v = {2tt)^ 1 J d 2 £yfgg ab (9 t Ml/ ) {dbt^ v ) is the number of self-intersections of the world 
sheet S. Here we have adapted the standard notations: g ab — {d a x^(S)){d b x lJL {$)) denotes the 
induced metric tensor of the surface, g = det ||g afc ||, and = e ab (<9 a £ M (£)) (dbX„(t;)) j \fg 
stands for its extrinsic curvature tensor. The coupling constant c can be calculated analogously 
to refs. |l^,^5) and reads c = — y|g (^-) ■ Since the obtained string 6*-term is proportional 
to the number of self-intersections of the world sheet, it might be relevant to the solution 
of the problem of crumpling of the world sheet, as it has been first mentioned in ref. p4| . 

2 

In particular, in the extreme strong coupling limit g — ► oo, c — > — yjg, and therefore self- 
intersections are weighted in the partition function with the desired factor (—1)" at 9 — ir/12. 

Notice also that owing to the field-theoretical 0-term, in this limit the dual mass ( |l5| ) is 
very large: 



Sir 2 ' 



(16) 



This situation is opposite to the case without #-term, where the dual mass equal to rjj g vanishes 
in the extreme strong coupling limit, since 77(g) oc e _const / s " — > 1. Therefore, contrary to that 
case, at 9 7^ 0, the expansion of the ^-dependent part of the action ( |i~4| ) in powers of the 
derivatives w.r.t. £ a 's [|13|, [L5[|, being nothing else but the expansion in the inverse powers 



( 5 ) The string 0-term can also be derived from the instanton gas model of QCD, as it has been done in 
ref. (li). 
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of to, converges fastly. In particular, this means that all the terms of this expansion higher 
in the derivatives than the Nambu-Goto, rigidity, and 0-ones are irrelevant, since they are 
suppressed by the higher powers of the mass. 

Clearly, all the lowest local string terms, i.e. the Nambu-Goto and rigidity ones, as well 
as the above-discussed 0-term, depend explicitly on the form of the world sheet S. This is 
the consequence of the fact that in a derivation of eq. (J13| ) we have taken into account only 
one (namely, real) branch of the monopole potential (in). However, as it has been argued 
first for the case of compact QED in rcf. || and then discussed for the 3D SU (3)-inspired 
case in ref . j^] , the world-sheet independence of the Wilson loop becomes restored upon the 
summation over all the complex-valued branches of the potential (^|) at every space-time point. 

Notice that although the low-energy limits of confining string theories obtained in compact 
QED and in the SU (3)-inspired case under study have the same functional form, the full 
expressions are different since in the S'[/(3)-case the B M „-field has two components. Moreover, 
the obtained results can be generalized to the S't7(-/V)-inspired case along with the lines of 
ref. ||. There the fields and have (N — 1) components, but the functional form of 
the action remains the same modulo the fact that the fields j( Q )'s are given by more 

complicated combinations of the fields |j^| a 's. Let us finally discuss the asymptotic behaviour 
m(N) at N ^S> 1. Note first of all that the strong coupling limit under consideration is still 
accessible at large N despite the fact that g ~ 1/ \/N . That is because the strong coupling 
limit implies only that g should be larger than some critical value, which itself scales as 1/y/N. 
Taking into account that r/(g) oc e - const /9 2 ; we see that according to eq. (15), the dual mass 



behaves in the large- TV limit as ViVe _const N . Besides that, in the SU(N)-case there exist 
N(N — l)/2 vectors q a 's with a > and consequently the same amount of mass terms of 
(N — 1) fields, which emerge from the expansion of cosines in the respective sine-Gordon 
theory describing these fields. Therefore, the number of mass terms of every field is of the 
order of N, i.e. each mass scales as \ N owing to this fact. These two observations considered 
together lead to the conclusion that at N ^> 1, m ~ N e~ const N . 

However, in the large-iV limit, the model under study becomes less and less relevant 
to QCD. That is because the number of off-diagonal fields disregarded within the Abelian 
dominance hypothesis, equal to [N 2 — N) , significantly exceeds in this limit the number of 
diagonal fields kept, equal to (N — 1). The problem of accounting for off-diagonal degrees of 
freedom deserves special investigations and will be considered in future publications. 



* * * 
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Appendix A. Details of integration over the Ti^-field in eq. (f7j[h - The desired integra- 
tion is Gaussian and can obviously be performed by finding the inverse to the quadratic part of 
the B Mt/ -ficld action. Since the integral under study factorizes into the product of the integrals 
over .B^'S) we should only solve the following equation in the momentum representation: 

1 fp 2 . 1. iS \ 
Here, the following projection operators were introduced (see e.g. ]lt|): 
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where "P M „ — <5 M „ — p^p v /p 2 ■ In what follows, we shall for brevity omit the indices and use 
the notation (O ■ O') x = O liVa j30' a p Xp . Then, representing G as G^ 1 ' + iG^ 2 \ we get the 
following equations for the real and imaginary parts of eq. (A.l): 



(V±P + . G (D + ° . G ( 2 ) = 21, ( P - 2 P + liV G^ - * . G M = 0. 

Applying further to the second of these equations the operation e-, it is possible to resolve 
it w.r.t. G^\ Substituting the so-obtained expression into the first equation, we get the 
following equation for G^ : 



■3- r 



P-l 



(P-i) 



G (2 



■G^ 



-1. 



Finally, seeking for G^ in the form of the Ansatz G^ 



P - 1 



e + e-(P-i\ 
f 



64tt 4 " 2tt 2 
f(p)s and using the equation 



e = 

2^2 



-41, we find 



<V?7 



G (1 



p2 _|_ m 2 



[?('-*) 


1 






V 1 



8tt 2 (p 2 + to 2 ) ' 

where the mass m of the field B^„ is given by eq. (|l5|). One can now straightforwardly obtain 
eq. ( |i"4| ) of the main text by passing back to the coordinate representation, using Stokes 
theorem, and the fact that for every a, Q 2 =1/3. 
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